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Abstract 

The estimates of overall strange quark production in high energy e + e~, pp and 
pp collisions by using the statistical-thermal model of hadronisation are presented 
and compared with previous works. The parametrization of strangeness suppression 
within the model is discussed. Interesting regularities emerge in the strange/non- 
strange produced quark ratio which turns out to be fairly constant in elementary 
collisions while it is twice as large in SPS heavy ion collision. 
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1 Introduction 



The interest for strangeness production in high energy collisions has been strengthened by the 
possibility of its use as a probe of the Quark Gluon Plasma formation [|IJ . An enhanced strange 
particle production in heavy ion collisions with respect to nucleon-ion collisions have actually 
been observed |2], |3| . In order to understand the relevance of such observations it is very im- 
portant to study strangeness production in high energy elementary collisions like e + e~or pp, as 
they provide a natural baseline to gauge the results in more complex systems. 
The comparison between strange and non-strange particle production in whatever kind of col- 
lision can hardly be accomplished by using only experimental data without any input from 
models. In fact, not all particle species yields are practically measurable; moreover, even if an 
experiment was able to measure all stable (disregarding weak decays) hadrons species, strong 
decays following the primary hadron production modify the number of u, d, and (to a less ex- 
tent) s valence quarks, thereby affecting the quantitative estimation of their actual production 
in the hadronisation process. The direct experimental measurement of strange to non-strange 
quark production therefore would require the determination of the inclusive yields of all hadron 
species, an exceedingly difficult experimental task. As a consequence, a model of hadronisation 
is necessary to estimate the production of a large number of unmeasured hadrons and reso- 
nances. Of course a prerequisite of any estimate is the agreement between the model and the 
existing data with a fairly small number of free parameters in order not to spoil its predictive 
power. 

In this paper I will discuss strangeness production by using a statistical-thermal model of hadro- 
nisation which turned out to be in good agreement with high energy collision data by using 
only three free parameters @,^|. 

2 Statistical model and strangeness suppression 

The statistical-thermodynamical model of hadronisation has been described and discussed in 
detail elsewhere [||,||. The main idea of the model is the existence of a set of pre-hadronic 
clusters (or fireballs) converting into hadrons according to the equiprobability of multi-hadronic 
phase space states (Gibbs postulate), where phase space is locally defined by the mass and the 
rest frame volume of each cluster. This assumption entails the possibility of expressing overall 
hadron multiplicities (i.e. measured in full phase space) by means of thermodynamical formulae 
within a canonical framework, that is with exact quantum number conservation. 
The model, whose basic parameters are the temperature T of the clusters/fireballs and the sum 
V of their proper volumes, is supplemented with an extra strangeness suppression factor, which 
is formally beyond a pure statistical hadronic phase space model. This factor 73 is introduced in 
the partition function used to derive hadron mulitplicities, as a fugacity related to the valence 
strange quark content of the hadron |7J . The canonical partition function supplemented with 
7s reads [|J: 

z(q0) = r d50 eiQ °' exp [v e s d3p iog (i ± 7^-v^/^)±i] . 

(i) 

where Q° = (Q, N, S, C, B) is a five-dimensional vector containing the initial electric charge 
Q, baryon number N, strangeness 5*, charm C and beauty B; Jj and rrij are the spin and 
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the mass of the j hadron species; q., is its quantum number vector (Qj, Nj, Sj, Cj, Bj) and Sj 
its number of valence strange+antistrange quarks. The hadron average multiplicities resulting 
from eq. (Jl|), in the Boltzmannian limit, reads: 

< n, > = -,' S >(2J, + dYL^K^) 2 ^-^ , (2) 

where the ratio Z(Q° — qj)/Z(Q°) is defined as chemical factor || and, unlike an ordinary 
fugacity, is not an intensive quantity for it depends on the system size. 

This way of regarding 73 is in fact a grand-canonical one, as fugacities can be defined only in 
a grand-canonical framework. On the other hand, it is possible to define 7s in a more general 
way which is independent of the adopted statistical formalism, i.e. microcanonical, canonical 
or grand-canonical. 

Indeed, for a certain multihadronic state {n\, . . . ,uk} (^1 is the number of hadrons belonging 
to species 1, is the number of hadrons belonging to species K) originating from the 

collision, its phase-space probability is multiplied by a factor 73 powered to the number of 
strange+antistrange quarks whose creation out of the vacuum is needed in order to set up that 
state. According to this definition, in a canonical framework, the probability P of realizing a 
multihadronic state {ni, . . . ,uk} in a collision whose initial state does not have any strange 
quarks is: 

P oc eM-E/T) T p K=1 " jSj 5^ . rij ^ Q0 , (3) 

where Sj is the number of valence strange+antistrange quarks contained in the j th hadron. This 
probability can be worked out to calculate the canonical partition function leading to the same 
expression as in eq. ([J). The above definition can be easily extendend to a microcanonical 
or a grand-canonical framework. Furthermore, it should be emphasized that this definition is 
more general and more appropriate for collisions with initial strange quarks (for instance K 
p) in which the use of the same partition function (TJ) obtained for colliding systems devoid 
of valence strange quarks would lead to odd results. In this special case, the probability (|3|) 
becomes: 

P oc eM-E/T) 7s E - 1 n ' sHSl 5 Ej njq3iQ0 , (4) 

where | S\ is the initial absolute strangeness, so that the canonical partition function will differ 
by a factor l/ 7 s 5 ' from the (P. 



3 Results on strangeness production 

The free parameters (T, V and 7 s) of the statistical-thermal model are determined by a fit to the 
measured hadron multiplicities (or multiplicity ratios) by taking into account the decay chain 
following the primary production. The quality of the fits is good and the multiplicity of up 
to 25-30 particles can be well reproduced. For heavy ion collision M one more free parameter 
is necessary (the baryon-chemical potential) as the number of participant nucleons, i.e. the 
total baryon number, is a measured quantity unlike in e + e~, pp and pp where it is known a 
priori. The extracted 73 values in elementary collisions at various centre of mass energy points 
are shown in fig. 1 along with those for heavy ion collisions at SPS energies M; also shown 
the result of a fit to AGS Si+Au hadronic ratios quoted in ref. (see table 1). All 7 s's turn 
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Table 1: Results of thermal model fit to particle ratios at AGS Si+Au collisions. The number of 
participant nucleons has been set to 108 and electric charge to 45 according to the estimation of 
ref. The fitted parameters are T, VT 3 exp[— 0.7GeV/T], 7s and Hb/T whilst X s is derived 
from the fitted parameters and their errors. The left column is the set of parameters obtained 
by setting the error on the number of participants to 10, the right column to 0.1 and keeping 
the same total electric charge. Besides experimental errors, also the uncertainties on masses, 
widhts and branching ratios of hadrons have been taken into account according to the procedure 
described in ref. |§ . The results are consistent with those in refs. ||, [U| . 



Fit 1 



Fit 2 



T (MeV) 

^T 3 exp[-0.7GeV/r] 


118.4±11.6 


132.8± 11.3 


0.73±0.25 


1.04±0.25 


7s 


0.923±0.16 


0.792±0.12 


X 2 /dof 


4.41±0.49 


3.89±0.34 


21.1/6 


22.5/6 



V (fm 3 ) 
As 



1250 



664 



0.488±0.093 0.569±0.10 



out to be definitely less than 1 except the AGS point which is consistent with 1, in agreement 
with the results of refs. |9|, [H]], though with a large uncertainty. This finding indicates that 
a strangeness chemical equilibrium at hadron level is not attained in most of the examined 
collisions. However, no apparent regularity emerges from the plot, as 7s has different values 
in e + e~and heavy ion collisions with respect to pp and pp collisions. The situation drastically 
changes when computing the ratio As: 



0.5((uu) + (dd)) ' v ; 

between newly produced valence ss pairs and half the sum of newly produced valence uu and 
dd pairs. As shown in fig. 2, As has a pretty constant value of around 0.2 for all examined e + e~, 
pp and pp collisions whereas it is twice as high in heavy ion collision both at SPS and AGS 
energies; the central As value at AGS is indeed larger than at SPS but the fit error is so large 
that a conclusion is not allowed. It should be emphasized that As is a strangeness suppression 
parameter related to the quark content of hadrons whereas 7s is related to the strange hadrons 
phase space. The computation of As proceeds as follows: firstly, the best fitted primary average 
multiplicities < rij > for all hadron species are used to count the total number of quarks < Q, t > 
of a given flavour i: 

<Q i >=^<n i >(T,V, Ts ) Qij (6) 

3 

where is the number of valence quarks of flavour i contained in the j th hadron. Secondly, 
the initial colliding quarks are subtracted so to count only the primarily newly produced ones. 
The errors on As are estimated by propagating the fit errors on T, V and 7s onto the primary 
multiplicities < nj >. Needless to say, the reliability of the estimate where only few hadron 
species were measured relies upon the agreement between model and data where a large number 
of measurements are available. 

The reason for the different behaviour of 7s and As in elementary collisions is twofold. On one 
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hand, a hadron gas with fixed 73, T and V has an increasing value of As for an increasing total 
baryon number (see fig. 3); hence, if As is constant, 73 must be lower in pp collisions (baryon 
number 2) than in e + e~and pp , provided that temperatures and volumes have similar values 
as it actually is ||. The physical reason of the increase of As with B in a hadron gas is the 
lower energy threshold for strange pair production in a baryon-rich environment where A + K 
production is favoured in comparison with kaon pair production in a baryon-free environment. 
On the other hand, if As is constant, a decrease of 73 with increasing V is expected, due to the 
so-called canonical strangeness enhancement HIT, 1~2] implying an increase of the strange hadron 
chemical factor (see eq. (|2|)) in increasingly large systems (see fig. 4). This effect accounts for 
the 7s difference between e + e~and pp collisions as the multiplicities and the volumes of the 
latter are somewhat higher for the considered centre of mass energy points. Other smaller con- 
tributions arise from the anticorrelation between 73 and T (outcoming temperatures in e + e~are 
slightly lower in e + e~than in pp ) and from the relatively more abundant production of heavy 
flavours in e + e~ collisions which slightly affect the fit of 7s. 

If As is more fundamental than 73, the best way of parametrising strangeness suppression would 
be to fix the mean absolute value of strangeness (with possible superimposed fluctuations) in- 
stead of using 7s . This would amount to extend the conservation laws from 5 (Q, N, S, C, B) 
to 6 (Q, N, S, C, B, \S\) in the calculation of canonical partition functions in eq. (|T|), where 
\S\ is now meant to be the overall number of strange quarks including both initial and newly 
produced ones. An interesting question, under current investigation, is why 73 has worked well 
anyhow in reproducing the strangeness suppression. This is not surprising for large systemes 
because the demand of exact \S\ leads to a fugacity, namely 7s (see for instance ref. fll3"l ), 
but it is a relevant question for elementary collisions where this is not possible; perhaps the 
replacement of 73 with As might lead to a further improvement of the fits. 
The As value estimated by the statistical-thermal model of hadronisation in elementary col- 
lisions is in agreement with a previous estimate quoted in ref. [[14] only for centre of mass 
energies y/s < 100 GeV, whereas its rise in pp collisions at \fs > 100 GeV claimed there is not 
observed. The reason of this discrepancy is the fact that As was estimated in ref. ]lj]| by using 
only K/tt ratio as experimental input and two parametrizations of hadron multiplicities |!5|,fL6[| 
which, unlike ours ||, do not satisfactorily reproduce all available measured multiplicities in 
pp collisions. As an example, for pp collisions at \fs = 546 GeV, the parametrization in ref. [|lj| 



predicts a A/K° ratio of 0.49, by taking the As = 0.28 value quoted in ref. [14j], whereas the 
experimental value is 0.24±0.05 



The enhancement of As in heavy ion collision at SPS energies is indeed very interesting. An 
explanation of this effect by advocating the formation of a fully equilibrated hadron gas is 
not viable, as 7s = 1 seems to be ruled out from recent analyses of the presently available 
data 0,[13|,|Ij|. On the other hand, it is not possible either to explain the rise of As by invoking 
a canonical strangeness enhancement from pp to heavy ion collisions at fixed 7s just because a 



fixed 7s would not be able to account for the production of meson in pp collisions |12 



A fully equilibrated hadron gas seems to be able to account for AGS data at \fs = 5.4 GeV, in 
agreement with the conclusions of refs. f9|,|I(]|, though the outcome error on 73 is quite large. 
However, two caveats for AGS fit should be mentioned; firstly, the ratios used in the fit have 
been measured in limited rapidity intervals whereas the model calculations involve full phase 
space multiplicities. Secondly, the fit has been performed within a pure grand-canonical frame- 
work, while a canonical treatment should bring about a slight 73 rise because a strangeness 
suppression due to the finite volume is involved. However, for the large extracted V values (see 
table 1), this is expected to be a minor correction. 
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4 Conclusions 



A statistical-thermal model of hadronisation supplemented with a non-equilibrium strangeness 
suppression parameter 73 has been used to study primary strangeness production in elementary 
and heavy ion collisions at high centre of mass energy. The reliability of the model is based on 
its success in reproducing many particle multiplicities at several centre of mass energies, with 
only three free parameters. It is found that the ratio of newly produced strange to non strange 
quark As (see eq. (H)) is fairly constant for all elementary collisions. This constancy emerges as 
an universal feature of hadronisation along with the constancy of temperature || . Thus, while 
u and d quarks are produced in the same amount, strange quarks are produced in a fixed ratio 
with u, d, which is most likely related to their higher mass. All quark pairs generated before 
or during the hadronisation process eventually fill hadronic phase space according to a pure 
statistical law. Hence, the fact that 73 < 1 in a modified hadron gas framework, simply reflects 
the underlying more fundamental As — 0.2 parameter. In order to emphasize the role of the 
hadronisation model in extracting the As regularity, it should be noticed that pure experimental 
observables such as the K + /tt + ratio does not show up such constancy (see fig. 5). In fact, kaon 
yield is enhanced in e + e~ collisions with respect to the examined pp and pp collisions because 
of the contribution of heavy flavoured hadron decays. 

In heavy ion collisions at SPS energies, As turns out to be twice as large. This enhancement 
cannot be accounted for a hadron gas explanation and it is a subject of reflection. There have 
been some attempts to explain it in terms of Quark-Gluon Plasma formation [12[ 19fl. 
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Figure captions 



Figure 1 Fitted 7s in e + e , pp, pp and heavy ion collisions as a function of centre of mass energy 
(nucleon-nucleon centre of mass energy for heavy ions) . 

Figure 2 Estimated As in e + e~, pp, pp and heavy ion collisions as a function of centre of mass 
energy (nucleon-nucleon centre of mass energy for heavy ions). The two different estimates 
for pp collisions correspond to the two limiting cases of annihilation (lower) or survival 
(higher) of the 6 colliding valence quarks. 

Figure 3 The As parameter in a hadron gas at fixed T, V and 7s as a function of the number of 
initial protons (baryon number equal to electric charge). 

Figure 4 Chemical factor Z(0, 0, 1, 0, 0)/Z(0, 0, 0, 0, 0) of a neutral strange hadron (e.g. K°) in a 
completely neutral hadron gas as a function of the volume for 73 = 0.5 and three different 
values of temperature. The chemical factor approaches its grand-canonical limit, i.e. 1, 
for large volumes. 

Figure 5 K + /n + ratio in elementary collisions as a function of centre of mass energy. 
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